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ABSTRACT 

The formulae relating turbulence and other atmospheric parameters to back- 
scattered power for radar observations are reviewed. The paper considers 
primarily the case of scatter from turbulent irregularities which have scales 
corresponding to the range of isotropic, inertial range turbulence, although 
some brief discussion of the applicability of this assumption is given. A new 
formula is introduced for the mesosphere which relates ionospheric electron 
densities to backscattered power. 

INTRODUCTION 

Discussions and the relationship between the intensity of turbulence and 
backscattered radar signal strengths have, in recent literature, been largely 
based upon the Kolmogorof f theory of inertial range isotropic turbulence (e.g. , 
BATCHELOR, 1953; TATARSKI, 1961, 1971). This is not to say, however, that this 
is the only possible approach. For example, BOOKER and GORDON (1950) and S TARAS 
(1952) adopted an alternative procedure for examination of turbulence (e.g., see 
review by GAGE and BALSLEY, 1980 ) . This second approach has not been as 
extensively applied as that due to Kolmogorof f, but, as pointed out by GAGE and 
BALSLEY (1980), it does allow extensions to conditions of anisotropic turbu- 
lence. Whether in fact the assumptions of inertial range, isotropic turbulence 
are valid for the atmosphere is to some extent an unresolved topic. For example 
the inertial range theory requires high Reynolds numbers (BATCHELOR, 1953, 
p 116), and Reynolds numbers in the atmosphere tend to be only moderate. 
Furthermore, observations of turbulence in the stratosphere often show very thin 
( *50-200 m thick) well-defined layers of turbulence (e.g., CRANE, 1980). This 
is not predicted by the Kolmogorof f theory. Nevertheless, it is normally 
assumed that Kolmogorof f theory still applies within the layer, at scales 
smaller than the layer thickness. BOLGIANO (1968) has proposed a turbulence 
model in which thin well-mixed layers of turbulence form, and in this model 
radio-wave backs cat ter is not produced by the turbulence within the layer but by 
discontinuities in refractive index at its top and bottom. The scatter from 
such discontinuities should be very different in character to turbulent scatter. 
It should show an aspect sensitivity, with most scatter coming from the verti- 
cal, and should have slow fading times. Such "specular reflections" are well 
known to occur in the stratosphere at VHF (e.g., GAGE and GREEN, 1978; ROTTGER 
and LIU, 1978), but whether the mechanism proposed by Bolgiano explains these 
reflections has not been resolved. Other refinements to Kolmogorof f theory have 
been presented by some authors (e.g., HILL and CLIFFORD, 1978; WEINST0CK, 

1978a). 

Despite these potential problems, however, the Kolmogoroff theory of 
inertial range turbulence appears to model the atmosphere reasonably well in the 
range of scales for which it is applicable. Therefore this model will be the 
main one discussed in this short essay. 

A short introduction of the formulae of the inertial range theory will 
first be given, and then it will briefly be shown how these formulae extend to 
radio-wave scattering. Some discussion on the accuracy of these formulae will 
then follow. 
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It will be assumed initially that the radar looks at the atmosphere at an 
off-zenith angle, so that the role of specular reflectors can be ignored. The 
complexities introduced by specular scatter will not be discussed in detail; 
more complete discussions can be found in, for example, HARPER and GORDON 
(1980), and ROTTGER (1980a, b). 


INERTIAL RANGE TURBULENCE 

Atmospheric turbulence causes random fluctuations of various atmospheric 
parameters, such as density, velocity, refractive index, etc. The statistics of 
the turbulence is usually described using one of these parameters. However, the 
parameter chosen to describe the fluctuations must be a passive tracer. This 
means that its statistical properties must not depend on the position in the 
turbulence patch. For example, density is not a good passive tracer, as dis- 
placement of a parcel of air vertically alters its density. This matter was 
discussed more deeply by TATARSKI (1961), and will also be considered again 
shortly. Potential temperature is a good tracer. So are the velocity fluc- 
tuations. 


For the present, let this passive tracer by a scalar, denoted by 9. 


There are at least two ways to describe the statistical properties of the 
turbulence. One way is by means of the structure function, viz 


D 0 (r) = < l e (x) - 9(x+r)| 2> , 


( 1 ) 


Here, x represent a position vector, and r. a spatial displacement. <> repre- 
sents an average over space and time. It can be shown that for inertial range 
turbulence. 


Vi) = c e 2r2/3 (2) 

2 

e.g. TATARSKI (1961), where C Q depends on the intensity of turbulence. The 
turbulence fluctuations can also be expressed as the Fourier sum of wave numbers 
k = 2 tt/A, A being the Fourier scale. Then TATARSKI (1961) showed that the 
spectrum of fluctuations is 

<f> 0 (k) = 0.033 C 0 2 |k| “ 11/3 (3) 

» « 00 

A normalization has been chosen such that JJf <I> Q (k)dk = <0 2 >. 

“ 2 

It can be shown (TATARSKI, 1961) that Cq is related to the outer scale of 
turbulence, L Q (i.e. , the approximate transition scale between the inertial and 
buoyancy ranges of turbulence) by the relation 


C 2 = a 2 a’L 

0 i 


4/3 


<f > 2 


(4) 


Here, a is a constant, = 2.8 (e.g^_, VANZANDT et a!., 1978), a* is a constant_ 
which is approximately 1., and (d<f/.dz.) is the gradient of the mean quantity 0. 


The formulae (3) and (4) form the basis of theories which relate back- 
scattered radar power to turbulence. However, before discussing how this is 
done, some other spectral forms should be briefly discussed. It is important to 
note that the spectral form shown in (3) is not the only form which appears in 
the literature. It is the full three-dimensional spectrum. But at times the 
spectrum of wave numbers with magnitude k *= |k| is given viz 

E 0 (k) = 4irk 2 (j> e (k) = 0.132it C 0 k~ 5/3 . 


(5) 
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No distinction between the directions of the k vectors is made in this formula. 
Another very important spectrum is the spectrum of fluctuations which would be 
seen by a probe moving in a straight line through the turbulence. This is not 
the same as (3), since that only looks at 1 scale direction. But for a probe, 
all scales produce an effect along the path of the probe, but their "effective 
scales'* change because they are not all orientated along the probe path. Then 
this produces a spectrum 
00 00 

s 0 (k) = / / * 0 (k)dk dk z (6) 

—00 — oo * 

whence S fk) =0.25 C 2 k -5/3 (7) 

8 8 

OTTERSTEN (1969) has emphasized the difference between (3) and (7), and pointed 
out that (3) is applicable for radar experiments, whilst (7) is applicable for 
in situ measurements. The two expressions should not be confused. 


POTENTIAL REFRACTIVE INDEX GRADIENT 


In considering backscatter of radio waves from the atmosphere, it is of 
course necessary to look at the fluctuations in refractive index caused by the 
turbulence. The refractive index n of air at centimetre and metre wavelengths 
is (TATARSKI, 1961). 

n = 1 + 10~ 6 x (79/T) x (p + 4800 e/T) (8) 

where T is temperature (absolute), p is pressure (millibars) and e is the water 
vapour pressure. However, it is more convenient to express n as a function of 
the potential temperature H and the specific humidity, q [= e/(1.62p)], Then 

n - n(p, H, q) (9) 


But n here is not a good passive tracer. This can be seen by the following 
consideration. Suppose an eddy moves from height to a height z 2 . Suppose 
that at the eddy was in equilibrium with its environment, and that at this 
height there was pressure pj, potential temperature H and specific humidity q T - 
Suppose that at height Z 2 , the atmospheric pressure i8rp 2 , and the environmental 
H, q are H 2 and q 2 . However, at z 2 the parcel itself has H = H^ and q = q^, 
since it is assumed to have moved adiabatically . Of course the pressure in the 
parcel is now P 2 * Hence the difference in refractive index between the parcel 
and its environment at z 2 is 

An = n(z 2 , p 2 , H 1> - n(z 2> p 2> H 2> q 2 ) 

parcel environment 


A ,9a 8H 3n . 
An = + 3? AZ 


( 10 ) 


where Az = - z^. 

This is not simply the difference in refractive index at heights z and z , 
which would be ^ 


An 


riHlH + 

dp dZ 


. 3n 3H , 3n 3q x . 

+ Az 


(ID 


The formula discussed in the previous section can be applied for refractive index 
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with 0 = n, but in equation (4), the term (d0/dz) should not simply be (dn/dz) 
as given by (11), but rather, from (10), 

dn __ x 3n 3H 3n 3q ^ /■-» o\ 

<dz ~ ^3H 3z 3q 3z' 1 ' 

This quantity is often denoted by M, and is called the generalized potential 
refractive index gradient. 

Fdr metre and centimetre scatter from the un-ionized atmosphere, (TATARSKI, 
1961), 


= -79xl0~ 6 p 
T 2 


• (1 + 


15,500 
T ; 


(f + r 

dz a 


7800 dcu 

d-A^OO) * dz) 


(13) 


The term r is the adiabatic lapse rate, 
a 

In the stratosphere and mesosphere, q = 0. However, once heights of greater 
than 50-60 km are reached, scatter from turbulence is enhanced by the existence 
of free electrons (ionization), and in this case M needs modification. 


For these circumstances 


n — n(N , Vjg) (14) 

Where N is the electron density, and v m is the collision frequency of electrons 
with neutral particles. (Pressure and temperature fluctuations also produce 
weak changes in n, as for the troposphere and stratosphere, but these effects 
are very weak compared to the effects of N and v^, and so can be ignored). 

HOCKING (1980, 1981) has shown that the appropriate generalized refractive 
index gradient for the ionospheric D region is given approximately by 



[I ( <£ + 

l T Mz 


V 


dN , N d£ ' , 3n_ _^m 
dz p dz 3v 3z * 

r m 


where p is the neutral air density. 
For the region 0-120 km, this equals 


M 


e 


Jn| f dT + n _|[ 
3N L T Mz a' dz 


(1.4x10 4 )N ] + 


.. 3v 
3n m 

3v 3z 
m 


At VHF in the D region n is related quite simply to N by the relation 


2 i - 1 » , 2 

n =1 -it r NX 
e 


(15) 


(16) 


where r £ is the classical electron radius, and X is the radar wavelength, 

3n/3v - 0, and 

m 


3_n 

3N 



X 


2 


so 


At HF and MF, the relation between n and N is more complex (e.g., BUDDEN, 1965). 
Thus M depends on the potential temperature gradient, the electron density 
gradient, and the neutral atmospheric density gradient. 
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In deriving (15), it was assumed that when a parcel of ionosphere is 
displaced, the ratio of electron density to neutral density remains constant, 
and that no change in the photochemical reaction rates occurs during such a 
displacement. HILL and BOWHILL (1979) have suggested that this might not be 
exactly true, but nevertheless equations (15) and (16) should provide a rea- 
sonable estimate of M . 

e 

SCALES OF THE INERTIAL RANGE 

Before proceeding to show how these turbulence formulae relate to radar 
backscatter, it is important to illustrate over what scales they can be - 
applied. 

At very small scales, the kinetic energy density contained by the eddies 
is diminished due to viscous effects, and much of the turbulent energy is dis- 
sipated as heat. This small scale range is often called the Viscous range”. 

At very large scales, buoyancy effects become important, and turbulent eddies 
taken on a ”pan-cake”-like appearance, with horizontal scales much larger than 
their vertical dimensions. 

An important scale for determining the boundary of the inertial-range to 
viscous range transition is the Kolmogoroff microscale, defined by 

n = (v 3 /e) 1/4 . 9(17) 

Here, v is the kinematic viscosity, and e is the turbulent energy dissipation 
rate. This is a scale well within the viscous range. The scale 

i = 7.4n (18) 

o 

is known as the "inner scale" (e.g. , HILL and CLIFFORD, 1978) and defines the 
approximate transition scales between the inertial and viscous ranges. (The 
constant 7.4 is only relevant for air.) 

The scale for determining the transition region between the inertial and 
buoyancy ranges is (WEIN STOCK, 1978b) 

L b = (2tt/ 0.62) e 1/2 w B _3/2 , (19) 

where is the Brunt-Vaisala period of the atmosphere at the height of the 
turbulence. This 1^ should not be confused with L in (4): they are different 
parameters, as will 13 be seen later. 0 

The inertial range of turbulence strictly only applies for scales somewhat 

less than L^ and larger than l . 

B o 

Approximate values of L g and are shown in Figure 1. For a radar wave- 
length \ backscatter occurs for scales of V 2. Thus if a 50 MHz radar is used 
(X = 6 m) , then scatter should be possible from isotropic inertial range 
turbulence up to altitudes of about 65-70 km. Figures similar to Figure 1 have 
appeared elsewhere in the literature (e.g., GAGE and BALSLEY, 1980 ), and show 
similar values for n and l * 


RADIO-WAVE BACKSCATTER 

Having illustrated some appropriate formulae for relating refractive index 
fluctuations to turbulence parameters, it is now necessary to determine how 
these formulae relate to radio-wave backscatter. 
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Figure 1. Typical inner and outer scales (& 0 and Lg respectively) for inertial- 
range turbulence in the atmosphere. The formulae used are given in the text. 

The profile of the Brunt-Vaisala period is also shown. It was assumed that 
the mean value of e was 10"^ W kg~l at 90 km, decreasing exponentially to 10~2 
W kg~l at 80 km. Between 80 and 60 km, this mean was taken at 10~2 w kg“l. 

In the region 60-90 km, the bounds of the dotted areas correspond to turbulent 
energy dissipation rates of 1/3 rd and 3 times these mean values. Below 40 km, 
and downto the tropopause, the upper and lower limits of e were taken as 10“ 3 
and 10~5 w kg~l. Lg and were assumed to vary smoothly between 40 and 60 km. 
(The region between 30 and 80 km is perhaps the most uncertain part of the 
graph.) Below the tropopause, e was assumed to be limited between 10“4 and 
10" 1 W kg~l. Larger e values correspond to smaller l Q and larger Lg values 
(i.e. the inertial range widen at both ends, when e increases). 


BOOKER (1956) has shown that the power backs cat ter ed from refractive index 
fluctuations with mean square value N(k)dk at scales A - 2n/k, per unit solid 
angle, per unit incident power density, and per unit volume (i.e., the cross 
section of backs cat ter) is 

a= [4 TV 2 /A 2 ]N(k). (20) 

(this expression is true at VHF, but at lower frequencies may not be) P(k) is 
similar to 4> n (k) in equation (.3) (0-n) , but BOOKER (1955) used the normaliza- 
tion 

(2tt) — 3 /// N (k) dk 1 dk 2 dk 3 = | An | 2 
—00 

Hence N(k) = (2 tt) 3 4> (k). (21) 

— n • — 

For radar backs cat ter at wavelength X, k « 4^r A and so using (21) and 
(29), and using (3) for <f> (k), 

, - .00654, 4/3 C 2 r 1/3 

n 


(22) 
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Sometimes an alternative definition of backscatter cross section is used. 
This is the total power which would be scattered if power were scattered iso- 
tropically with an intensity equal to that of the backscattered radiation, per 
per unit incident power density, per unit volume of scatterer. It is often de- 
noted by rj, and although this can at times be confused with the Kolmogoroff 
microscale, this convention will be maintained here. Then 

q = 4 tt a, (23) 


and hence 


n = 0.38C X 
n 


*1/3 


(24) 


as also derived by OTTERSTEN (1969), It should be noted that, at least for the 
ionosphere, the wavelength dependence for q is more complex than X~l/3, since 
C^ 2 is itself a function of wavelength. This can be seen by considering a 
given patch of turbulence in the ionosphere. For this patch, there is a con- ^ 
stant electron density structure constant, which we may denote by C Then C 
= (an/3N) 2 C 1 (e.g., HOCKING and VINCENT, 1982), and (3n/8N) is strongly N 
wavelength dependent, as has already been seen. 


Now it is necessary to show how a (or q) relates to the power received in a 
backscatter experiment. Consider scatter from a height h. Then the peak power 
per unit area incident at h is 


p h = (P T G T e T )/(4,h ), 


(25) 


where is the transmitter peak power, G^is the transmitter array directivity 
and e T is the transmitter efficiency. To obtain the power backscattered per 
unit steradian, we simply multiply this by aV, where V (the "radar volume") is 
the volume defined by the locus of the half -power points of the radar, and the 
pulse length. The receiving array subtends a solid angle A R /h 2 to this scat- 
tering region at height h, so the peak power received by the receiver is 
(ignoring absorption) 


P 


R 


P T S Vg e T P R A R 
4tr 


(26) 


where e R is the receiver efficiency and A R is the effective area of the 
r e ce iving a rr ay , 


Thus by (22) and (26), 

„ 16tt 2 P r h 4 X 

C 2 = 

n 0.38V P G e 


1/3 

e R \ 


(27) 


2 

In the case of a circularly symmetric array, V ~ iKhQ-j/o) • 
pulse- length and Qt_/2 ** s t ^ ie k a lf~P° wer ^h a lf“ w i-dth -Ot the polar 
the same array is used for transmission and reception, e T 55 e p , 
4ir, so 


where L is the 
diagram. If 2 
ahd A r = G t X / 



n 


2 -5/3 

1662.25 P R h ' X 1 
P T S' e T 2 e i/2 2 L 


(28) 


2 2 

For a radar, it can be shown that G T 6^ 35 tt /4, where 0^ is the half-power- 

ha If -width. For the case in which the same radar is used for transmission and 
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reception, |he effective half-power width reduces by / 2 , so /2 =0 l/ y ^» an< * 

G T 9 l/2 2 = IT /8. Then 

128 A 1/3 h 2 P 

C 2 = -5 2- (29) 

0. 38tt e A IP 

2 

Either (30) or (31) can be used to estimate C n from absolute measurements of 
received power. This has been done, for example, by NASTROM et al. (1982) and 
GOOD et al. (1982). 


APPLICATION OF THE FORMULAE IN THE REAL ATMOSPHERE 

2 

Clearly (28) or (29) can be used to estimate C n with a radar, but 
unfortunately this does not give the structure constant C n z for the turbulence 
itself. It would, if the turbulence filled the radar volume, but in reality 
turbulence appears to occur in thin horizontal layers, with depths of 10 s to 
100s of meters (e.g., VANZANDT et al., 1978; CRANE, 1980; WEINSTOCK, 1981, and 
references therein). Thus the scattering within the radar volume is usually 
from a few thin turbulent layers, and the effective volume V should not be 
1 T h 2 ®i/ 2 2 L as proposed earlier, but ('irh^G-j^i’) where F represents the 
fraction of volume within the radar volume which is filled with turbulence. 
VANZANDT et al. (1978) obtained a formula enabling F to be determined from a 
knowledge of the mean wind shear (taken with a resolution of about a kilometer 
or so), the* standard devi ation of the fine-scale shear, and a "critical wind 
shear” S c . Then, if C ^ 2 (turb) is the refractive index structure constant for 
the turbulence, and C n 2( radar) is the value measured by the radar, 

C n 2 (turb) * C n 2 (radar) /F . (30) 

VANZANDT et al.(1978) used meteorological data to estimate M (equation 13) 
and then applied (4) and (30) to estimate C n 2 ( radar). They found that with a 
value of Lq equal to 10 m, good agreement occurred between the model and radar 
observations, particularly in the stratosphere. Agreement was not so good in 
the lower troposphere, because the model did not account for humidity fluctua- 
tions. 

VANZANDT et al. (1981) improved the theory of VANZANDT et al. (1978) by 
considering these humidity fluctuations, by considering small-scale fluctuation 
in potential temperature, by letting the layer thicknesses be non-constant, and 
also by using* a more realistic distribution for the wind shears. 

* ___ 

; 2 

VANZANDT et al. (1978) compared C n (radar) estimates from their model to 
experimental radar values, assuming Lq 5 * 10 m. GAGE-et al, (1980) applied this 
principle' in reverse, using radar estimates of G n 2 to effectively estimate L 0 
(through equation 14). They then drew on an equation relating L 0 and the 
turbulent energy dissipation rate e, to estimate e. This relation was 
(TATARSKI 1961) 

. s = b's 3/2 L 2 (31) 

■ ° 

- 2 

where b was tafen as a constant equal to unity, and S - (du/dz) is the shear 
in |he mean wind. By replacing S. with o) 2 B/R^, where is the Richardson 
number and ujg is the Brunt-Vaisala frequency, and assuming that turbulence 
exists if = R Xcrit ), (a critical value), they obtained, using L 0 from (4), 

e , = [C 2 (turb) • (a^a VR. (crit) a>- 2 M 2 ) ^] 2 ^ 2 . (32) 

turb n l B J 

They took R i( cr ^ t ) * i/4. Then e tur jis the mean turbulent energy dissipation 
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rate, GAGE et al. (1980) also calculated a quantity which they denoted by e", 
which was the mean turbulent energy dissipation rate averaged over the radar 
volume* They took 


T = F 1 e 


turb 


(33) 


2 

From radar measurements GAGE et al. (1980) estimated ,C n (radar). Then 
they made some__reasonable assumptions concerning F, and so were able to esti- 
mate £ turb ^4 e from their radar data. M was calculated from meteorological 
measurements of T and p, and it was assumed that the humidity terms in M were 
unimportant. It should also be noted that F is dependent on cog , although this 
may not be obvious in the simplified discussion given above. This dependence 
of F on can cause some problems in estimating F, but GAGE et al. (1980) were 
careful to reduce this error as much as was reasonable. 


The technique outlined above is, at least in principle, the primary means 
by which e is obtained for the atmosphere using VHF radars. Variations on the 
details of these formulae have been presented (e.g., CRANE, 1980; WEIN STOCK, 
1981), but the principle remains similar - namely, to determine the fraction of 
the radar volume actually filled by turbulence, and then to correct C n ^ values 
measured by the radar to give e turt> and e. 

In their model calculations, VANZANBT et al. (1978) chose L 0 ^lO m. It 
should be noted that h Q is not equivalent to Lg in (19). If S = Wg /R is 
substituted in (31), as was proposed, then 

1 Q = 0,35 e 1/2 Ub -3/2 ( 34 ) 

and comparison with (19) shows that 

L 0 « .035 L b (35) 

The difference arises because of the different definition used to define 
these "outer scales "* Lg is probably a better measure of the transition scale 
between the inertial and buoyancy subranges, but L 0 is quite appropriate 
wherever the formulae of TATARSKI (1961) are applied. This of course means 
that (31) and (4) are only applicable for L Q as defined by TATARSKI (1961). 

The choice of L 0 « 10 m used by VANZANDT et al. (1978) corresponds to a choice 
of Eg of about 290 m. WEINSTOCK (1981) developed his theory relating £ and 
C n 2 ( radar) using Lg as an estimate of the sum of the thicknesses of the 
turbulent layers in the radar volume, and achieved numerical results similar 
(to within a factor of 2) to those of GAGE et al. (1980). 


It is also possible to apply (32) for the mesosphere, using M e (equation 
16) in place of M. However, there are some problems in estimating F for this 
case. For example, CZECHOWSKY et al. (1979), using a 150 m resolution radar, 
have shown that at mesospheric altitudes of »80 km, the scattering layers can 
be quite thick (up to 3 1 km) and so F may approach unity. ' Further, at VHF the 
appropriate scattering scales may be within the viscous range, so (3) and hence 
(29) may not be applicable. At HF and MF radar wavelengths (e.g., X = 150 m) , 

however, scatter should be from the inertial range and these formulae should be 
appropriate. 


DISCUSSION 

Interestingly, RASTOGI and BOWHILL (1976) presented some formulae relating 
turbulence parameters to backscattered power, and concluded that for the meso- 
sphere the backscattered power was independent of e. They based their con- 
clusions on dimensional arguments* However, these arguments were nowhere near 
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as rigorous as those presented in this paper, and it is felt that (32) more 
appropriately represents the relation between e and C n 2 . The appropriate 
generalized refractive index gradient can be obtained from (13) (troposphere 
and stratosphere) or (16) (mesosphere), provided that scatter can be assumed to 
be from within the inertial range of turbulence. It will be noted that (32) is 
indeed dimensionally correct, since C n 2 has units of m~ 2 ' 

Direct, independent measurements of e and (^(radar) have not been exten- 
sively made, so it is difficult to comment on the validity of these theories. 
Certainly, however, the estimates of e presented by GAGE et at. (1980) are of 
the correct order of magnitude. 

Recently, HOCKING ( 1983a, b) has presented an alternative method for 
measurement of turbulent energy dissipation rates with radars. This utilizes 
not the signal strength backs cattered, but the spectral widths of the received 
signal. The principle of the method has been known for many years (e.g., 

ATLAS, 1964; FRISCH and CLIFFORD, 1974; FRISCH and STRAUCH, 1976) but the major 
advance presented by HOCKING ( 1983a, b) was the accurate removal of both (i) 
spectral broadening due to the motion of the mean wind across the finite beam- 
width and (ii) spectral "broadening” (or "narrowing” in some cases) due to 
vertical wind- shears in the horizontal wind. These two factors have previously 
been considered separately (e.g. , ATLAS, 1964), but never coherently, HOCKING 
(1983a) also illustrated that there was a necessity to distinguish between 
vertical and horizontal fluctuating motions, and showed how this could be done. 
This technique was illustrated using an BF radar to measure energy dissipation 
rates in the vieso sphere. 

The technique can readily be applied at VHF, and the author is currently 
doing this with the "SOUSY" radar (ROTTGER et al., 1978) in West Germany. The 
estimates of e appear to be of the correct order of magnitude, and will be 
reported separately in a later paper. Application of this method can allow 
independent comparisons of e and C n 2 (radar) , and therefore can check the equa- 
tion (32). This new method of obtaining involves less assumptions than (32), 
and may prove to be a more reliable method in the future. Previously, some 
authors made comparisons of signal fading time and received power (e.g., FUKAO 
et al., 1980a, b; RASTOGI and BOWHILL, 1976b), but it is difficult to decide 
how much the fading time (or equivalently the spectral width) is contaminated 
be beam- and wind-shear broadening. Therefore these measurements cannot really 
be regarded as comparison of e and C n 2. 

For approximate estimates of the effects of beam-broadening, the following 
formula may be useful. If 8^/2*- 8 the half-power-half-width of the effective 
radar beam, and V is the mean velocity of the scatterers tangential to the beam 
(usually this amounts to the horizontal velocity), then the half-power spectral 
half-width due to beam broadening is 

f 1/2 - d,o) 2/x e 1/2 v (36) 

This is very nearly exact, provided beam widths of less than 3° - 4° are 
used. A similar formula was presented by ATLAS (1964), and was originally 
derived by HITSCHFELD and DENNIS (1956). [Atlas gives an equation ct v - 0.3 0V. 
However, this equation is for the case in which the same radar is used for 
transmission and reception, and 8 is the half-power- full-width for the trans- 
mitter (or receiving) polar diagram only. Thus 8 in this equation is equal to 
2/2 times B^^in (36)., since S^^there is the half-width for the effective 
polar diagram (transmitter and receiver polar diagram included).] For proper 
removal of beam-broadening and shear-broadening, however, the complete treat- 
ment presented by HOCKING (1983a) is recommended. 

As discussed in the introduction, there may be problems with the as sump- 
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tion of inertial range turbulence, and it is useful to list some of these. 
Specular reflection has already been mentioned, and the cause of this has not 
been fully explained. It is sometimes assumed to be a process separate from 
turbulence, but this may not be. For example, the model proposed by BOLGIANO 
(1968), which was discussed earlier, may be important. In this case, tilting 
the radar beam from the vertical may cause the layer to disappear, since one of 
the assumptions of Bolgiano's model was that turbulence mixes the layer so well 
that no parameters such as density vary with height within the layer. Thus the 
generalized refractive index gradient within the layer is close to zero, and 
very little radio-wave scatter from the turbulence itself can occur. If such 
layers do exist, and are not seen by tilted VHF radars, this could lead to 
biases in estimates of ~s for the atmosphere. For vertically beamed radars, the 
relation between the specular scatter and the degree of turbulence may not be 
simple. Investigations of this matter await more experiments. The possibility 
that turbulence could be anisotropic even at scales of meters has also been 
briefly mentioned. Multi frequency experiments may help resolve some of these 
issues. 


CONCLUSIONS 

If it is assumed that radio-wave scatter is from inertial range turbu- 
lence, then the back-scattered power and the energy dissipation rate can be 
simply related through equations (28) (or 29), (30), (32), (16), (19) and (33). 
The derivation of these equations assumed inertial range isotropic turbulence, 
and the scales within which this is probably true are indicated in Figure 1. 
Considerable experimental work remains to be done to determine when these rela- 
tions are valid, and when they break down. 
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